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A Hamiltonian based approach using spatially localized projection operators is introduced to give 
precise meaning to the chemically intuitive idea of the electronic energy on a quantum subsystem. 
This definition facilitates the study of electronic energy transfer in arbitrarily coupled quantum 
systems. In particular, the decomposition scheme can be applied to molecular components that are 
strongly interacting (with significant orbital overlap) as well as to isolated fragments. The result 
leads to the proper electronic energy at all internuclear distances, including the case of separated 
fragments, and reduces to the well-known Forster and Dexter results in their respective limits. 
Numerical calculations of coherent energy and charge transfer dynamics in simple model systems 
are presented and the effect of collisionally induced decoherence is examined. 
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I. INTRODUCTION 



Measuring, predicting and ultimately controlling the rate of electronic energy transfer (EET) within and between 
molecules is currently of great scientific and technological interest. For example, technologically understanding natural 
EET processes [1-3] may well contribute to of our ability to design efficient photovoltaic devices [4-8]. 

Ideally, given a large molecular system or molecular network through which electronic energy flows, we would 
proceed as follows: After obtaining the system Hamiltonian H , we would solve the time dependent Schrodingcr 
equation to obtain the wavefunction ^(t) for a given initial wavefunction *(0). To track the electronic energy Ea 
on, e.g., site A, we would compute Ea — (\D r (£)|'Hyi|\Ij'(t)), where Ha is the operator which corresponds to electronic 
energy on site A. The computed Ea, being an integral over an operator on the electronic degrees of freedom, would 
automoatically include decoherence effects arising from other degrees of freedom [ ]. 

Clearly this is computationally difficult for realistic molecular systems. Nonetheless, one expects, from basic quan- 
tum principles, that there exists such a real valued operator Ha- For this reason it is surprising that we find no study 
of such an operator for electronic site energies. Rather, EET is most often studied via a series of approximate schemes 
reliant upon simplifying models that are applicable in restricted circumstances. 

In an effort to better quantify the rates of EET, the role of vibrational modes in decohering EET, the effects of 
assorted coupling scenarios in molecular systems, the role of orbital overlap between sites, etc. we obtain below a 
meaningful Ha operator that describes electronic energy on a site. The operator is valid at all internuclear distances 
and, as such, resolves problems (elucidated below) arising from electron antisymmetrization at all distances. In 
addition, the model is applicable in all electronic energy transfer regimes, i.e. weak, intermediate and strong electronic 
coupling, and reduces to known results, e.g. the Forster and Dexter results, in the appropriate limit. 

This paper is organized as follows: Section II motivates the problem of defining electronic energy on a site via 
the simple example of the Hydrogen molecule, H 2 . Section III then introduces the Hamiltonian decomposition that 
effectively resolves this issue. Limiting cases are discussed in Section IV. Numerical results of coherent energy and 
charge transfer dynamics between Hydrogen atoms and molecules are presented in Section V and a summary with 
conclusion in Section VI. 



II. MOTIVATION 

A central issue associated with defining an Ha arises due to electron antisymmetrization. For simplicity, consider 
a molecule composed of two atomic fragments A and B, where interest is in the electronic energy in Site A or B. 
Motivation arises from considering even the simplest of cases, i.e., hydrogen. Here, the Hamiltonian for HD (chosen 
to simplify the issue by using distinguishable nuclei) is of the form 

H = H A i + H B2 + V h2 (1) 

where the A and B label the two nuclei, and 1 and 2 label the electrons. The electronic Hamiltonian is (with h = 1) 

H = -iv? - — ~Vl ~— + ^ + - — ■ (2) 

2 r A \ 2 r B 2 R r 12 r A 2 r B i 



Hai Hb2 V\,i 

R is the internuclear distance, r a i is the distance of electron i from nucleus a, and Vf is the Laplacian associated 
with the i th electron. The exact wavefunction, 1*1,2), satisfies 

#1*1,2) =£1*1,2). (3) 

We highlight the definition problem by calculating the electronic energy Ea of subsystem A for an initially naive 
choice of the subsystem Hamiltonian at large A — B interatomic separation. Chemical intuition requires that the 
resultant Ea be equal to the isolated subsystem energy for any acceptable definition of Ha- We calculate this energy 
first by incorrectly assuming that the electrons in the system are distinguishable and then contrast the result with 
the calculation that takes proper account of electron indistinguishability through wavefunction antisymmetrization. 

A. Ignoring Electron Indistinguishability 

In the absence of antisymmetrization, the molecular wavefunction is, to a first approximation, a simple Hartree 
product 

\* 1>2 } = \A 1 B 2 ) = \A 1 )\B 2 ). (4) 
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where Hai\Ai) = E A \A X ) and H B2 \B 2 ) = E B \B 2 ), and \Ai) and \B 2 ) are normalized spin orbitals centered on the 
respective atoms. 

The energy of the molecule AB is then 

E= (*i, 2 \H\$ h2 ) =E A + E B + V AB , (5) 

where 

E A = (A X \H AX \A X ) = (A x \ - \v\\Ai) + (A x \ - — \A X ) =T A + V A 

2 r A i 

E B = (B 2 \H B2 \B 2 ) = (B 2 \ - \v\\B 2 ) + (B 2 \ - —\B 2 ) = T B + V B 

2 r B2 



V AB = {A X B 2 



AiB 2 ) (6) 



11 1 1 

R r 12 r A2 r B1 

As R — > 00, V AB — > since all terms are 0(1/ R). The results meet expectations: at infinite separation the energy 
associated with atoms A or B is equal to the energy of the respective isolated atoms. Hence it appears that Ti. A \ = H A i 
is a possible candidate for measuring electronic energy on site A. 

B. Antisymmetrized Electrons 

Consider now the case where proper account is taken of wavefunction antisymmetrization. Here, the wavefunction 

is 

\*t 2 ) = ~\A 1 B 2 -B x A 2 ) (7) 

where the subscripts denote the electron label on the indicated atom and the superscript "S" denotes properly 
antisymmetrized quantities. 

Once again, we naively define [Eq. (2)] the energy on atoms A and B,via Yi A \ and H B2 , and the interaction energy 
V AB , as: 

£iH*f 2 l#Ai|*f, 2 ) 

E S B = (*l >2 \H B2 \*t >2 ) 

Vl B = (*f, 2 |V ll9 |*f,>. (8) 



Expanding E^ gives 



E S A = ^[(AiB 2 \H A i\AiB 2 ) - (A 1 B 2 \H A i\BiA 2 ) 



which, in the large R limit, becomes 



- (B X A 2 \H A1 \A X B 2 ) + (B 1 A 2 \H A1 \B 1 A 2 )] . (9) 



E s A = l -{E A +T B )= l -{T A + T B + V A ) . (10) 



Similarly, 

1 , . „ „.,„.„ „ . , 1 



^1 = ^(^2-^21^21^152-^2) = ^(T A +T B + V B ) . (II) 



Finally, consider the interaction term 



VL = \{A X B 2 - B x A 2 \V h2 \A x B 2 - B X A 2 ) . (12) 

As R -> 00, V AB = V A + V B) where 

(B l A 2 \V ia \B l A 2 ) = f {B* X B X ) (- 1 —) {A* 2 A 2 )dn, 2 = V A + V B as R -> 00 (13) 

J \r X2 r A2 r BX J 
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Hence, as R — ► oo 

VL = \(Va + V b ) 
E S A = \(T A +T B + V A ) 

E S B = \(T A + T B + V B ) (14) 

Thus, for the case of the antisymmetrized wavefunctions, although the total energy 

E s =El + E s B + Vl B =E A + E B . (15) 

is correct, the energy partitioning amongst individual site energies E A and E B [Eq. (14)] is wrong. 

This issue is far from being new. For example, Margenau raised the same concern in considering antisymmetrization 
issues [10] and attributed its resolution to a sophisticated measurement problem. Alternative suggestions would be, 
for example, that antisymmetrization becomes unnecessary at large distances between electrons. This, however, is 
not the case. That is, requirements for electron exchange do not arise via a force that diminishes with R. Rather, 
examination of the proposed site operator 7^ shows, as discussed below, that the likely problem is that H A i,H B2 
themselves unjustifiably distinguish electron 1 from electron 2. 

III. DEFINING THE SUBSYSTEM HAMILTONIAN 

Consider now the general case. The electronic Hamiltonian of the composite system, H , may be written explicitly 
in terms of its one and two electron terms, h and g, as 



H 



N N 



where h{xi) — — ^Vf — Zp/rpi, g(xi,Xj) = l/Vy, i and j label the electrons, a labels the nuclei and Zp is the 
charge, in atomic units, on the j3 nucleus. The prime on the second summation denotes i =/= j. 

Since the electrons are indistinguishable and all interactions are two-body we can simplify the integrand when the 
Hamiltonian appears inside an integral over electronic coordinates as follows: 

(*(xi,x 2 , . . . ,x N )\H\$>(xi,x 2 , • • -,x N )) 

N N 

= (&(xi,X2, ■ ■ .,x N )\y]h(xj) + - 2J g[x i ,Xj)\'9(xx,X2, ■ ■ ■ ,%n)) 

i,j (17) 

N(N - 1) 

= (*(xi,x2, . . .,x N )\Nh(xi) H g(x 1 ,x 2 )\ 1 t>{x 1 ,x 2 , ■ ■ - ,x N )) 

= ($(xi,X 2 , . . .,X N )\H\^>(Xl,X2, ■ ■ .,x N )) 

where ^(xi, x 2 , ■ ■ ■ , xn) is a properly antisymmetrized wavefunction, and 

N(N - 1) 

U = Nh{ Xl ) + ^— 2 > -g{x l ,X2). (18) 

Giving Eq. (17) we can use Eq. (18) in place of the Hamiltonian whenever it appears inside an integral over all 
electronic coordinates. This substitution is also valid if the integrand is of the form HVL for any that is symmetric 
with respect to pairwise electron interchange and is composed of one and two electron operators. This is used below 
to simplify the notation. 

Without loss of generality we assume that our interest is in the electronic energy on a single subsystem, denoted 
A. The subsystems, each labeled by the index a, are defined spatially such that they collectively span all space. 
Given the system Hamiltonian H, we set out to define a subsystem Hamiltonian, 'H A , such that the electronic energy 
of the subsystem, E A , is given by E A = Tr[p% A ], where p = is the system density operator. In coordinate 

space ty(xi,x 2 , . . . ,xn) is the normalized antisymmetric wavefunction of the system, and the trace is taken over all 
electronic coordinates. Specifically, we look for a form Ha — ^aH, where T A is a superoperator. Interestingly, the 
choice of T A depends upon whether interest is in stationary or non-stationary states. 
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Three principles guide the choice of Ha and associated fragment energy Ea = Tt[pHa]'. 

(1) The fragment energy Ea must be real, for both stationary p as well as time dependent pit). 

(2) The energy Ea must reduce to the correct energy for an independent fragment A. 

(3) The operator Ha must be symmetric with respect to electron interchange. 

The first two requirements are evident, whereas the third benefits from some justification. Specifically, consider 

E A = Tv[pH A ] = (V(t)\H A \*(t)) (19) 

Since p is symmetric with respect to electron interchange, if Ha were antisymmetric with respect to electron inter- 
change it would result in an Ea which is similarly antisymmetric. Such an energy, whose sign would depend on 
electron identity, is non-physical. As a consequence, Ha must be symmetric. 

In addition, the method should reproduce standard results (e.g., Forster and Dexter energy transfer) when applied 
to those cases. All of these characteristics are verified for the Ha definition introduced below. 

Stationary state. The choice of V a in the definition of the subsystem Hamiltonian Ha is motivated by physical 
intuition in conjunction with the requirements above. In particular, the proposed expression for the subsystem 
Hamiltonian is 

N(N - 1) 

H A = TaH = r A [Nh( Xl ) + — ^ -g{x u x 2 )} 

N(N-l)^ ( 20 ) 
= NQ A (xi)h(xi) -I 2 J ®A,a(xi,X2)g(x 1 ,x 2 ) 



where 



Qa(xi) = { 1 . (21) 

otherwise . 



Here 6^,0(^1,^2) is a symmetric spatial projection operator 



8^(11,83) = \{Q A {xi)Q a {x2) + e a {x 1 )e A (x 2 )) (22) 



such that 



® a , A (xi,X 2 ) = &A,a(xi,X 2 ) 

9 A , a {xi,x 2 ) = Q a ,A(x 1 ,x 2 ). (23) 
Note the completeness of the projection operators: 

l(a;i) = 5^e a (ar 1 ) (24) 

a 

and that Li and Parr [11] have shown that &A,B(xi, Xj) is "physically sound" [12] and that it preserves the completeness 
relation which generalizes, for two electrons, to 

3:3) = 9g,p(xi,X2). (25) 

a, p 

Here a = A, B and j3 = A, B. With these definitions we can define Ha concisely as 

N(N — 1) N(N-l)x^ 

H A = Nh A {xi) H 9a,a{xi,x 2 ) H ^ 9A,ot(xi,x 2 ) (26) 

a^A 

where h A (xi) = Q A {xi)h(xi), gA,A(xi,x?) = &a,a(xi, x 2 )g(xi,x 2 ) and gA, a (xi,x 2 ) = A , a (xi,x 2 )g(xi,x 2 ). Hence 
hA{xi) is the one electron energy of the electron density that resides inside region A, gA,a( x ii x 2 ) is the two electron 
energy when only one of the (indistinguishable) electrons is in region A, and gA,A(xi, x 2 ) is the two electron energy 
with both electrons in region A. 

The resultant expression for the subsystem electronic energy, Ea, is valid for general time independent densities, 
p = Y^iPil^iti^ili giving the electronic energy on site A as 



G 



E A = Tr{H A p} = ^2piTr{U A \9i){9i\} = Tr 12 {H A p 12 } (27) 

i 

where p\ 2 = YliPi f dxsdx^ . . . da; J v(x|^i)(^i|x}, and x denotes (iei, a; 2 , . . . , xjy). The coordinate representa- 
tion of pi2 is (xi, x 2 \pi2\xi, x 2 ). In evaluating terms like Txx 2 {h A (xi)pi 2 }, we calculate Tr^/i^xi)/?]^} = 
I x ' 1 =x 1 dxidx 2 h A (xi)(x 1 , x 2 \pvi\x\, X2) , where we put 3^ = x\ and a; 2 = x 2 after operating with ^^(xi) but be- 

fore carrying out the integration [13]. That E A is real follows from the reality of the terms in the sum of products 
above. All the terms in the subsystem Hamiltonian T~i A (Eq. 26) are evidently real and, for a system in a stationary 
state, the coordinate representation of the density matrix is also real. 

Time Evolving State. For a system in a non-stationary state, significant for energy transfer studies, the wavefunction 
is described by a superposition of energy eigenstates. The matrix elements of the time dependent density operator, 
p{t) — ^iPi\^ i{t)) i{t)\, in the coordinate representation are generally complex. We therefore need to generalize 
r^4 to ensure that E A (t) is real. To examine the issue note that given a one particle complete orthonormal basis {xk} 
one can express the subsystem energy E A (t) as 

E A (t) = Tv{n A p(t)} = Tr 12 {HAPi 2 (i)} = J2 Pi™, jk nf lm (28) 

jklm 

where 

Pjk,lm = (x J (l)Xfe(2)|/5i 2 |x/(l)Xm(2)), (29) 

and 

nf lm = ( Xj (i)xk(2)\T A n\ X i(i)xm(2)) 

= N( Xj (l)\T A h(l)\ Xl (l))S km + N{N ~ 1} ( Xj (l)x fc (2)|r Ag (l, 2)lx,(l) Xro (2)) . (30) 

Consider each term in Eq. (28). Since p(t) is Hermitian, p* k im (t) = pi m ,jk(t)- While each term in the two electron 
contribution is generally complex, the fact that the projection operator, T A , and the two electron interaction operator, 
g(xi, X2) = I/X12 = l/\x\ — x 2 |, are real ensures that the terms in the expression for E A (t) sum in pairs to give a real 
result. In particular this is because ( Xj (l) Xk {2)\T A g(xx, X2)\ X i(l) Xm (2))* = { X j{l) X k(2)\T A g{xi,x 2 )\ X i(l) X m{2)) for 
any choice of region A subject to the conditions on the projection operators discussed above (Eq. 21-25). 

By contrast, the one electron term (xj(l)|r ' A h{xi)\xi{^)) contains a kinetic energy contribution T = — V^/2, 
which is Hermitian only for particular partitions of space such that the wavefunction satisfies particular boundary 
conditions. These boundary conditions are trivially satisfied when the wavefunction and its first derivatives vanish 
at the boundary. This is the case, for example, when we consider the entire system, where its boundaries are at 
infinity. The wavefunctions also vanish in the region between two subsystems if they are sufficiently separated. This 
is, however, generally not the case. 

To bypass this difficulty we recognize the possibility of utilizing non-Hermitian operators that have real eigenvalues 
[14]. We introduce an alternative Hamiltonian-based real space partitioning approach that is computationally efficient 
and applicable to arbitrary electronic states of the system. Specifically, we average the sum of the matrix with its 
transpose. The generalized subsystem Hamiltonian, Tl A , obtained in this way is applicable to both time dependent 
and time independent densities and is of the form: 

n, r n, AT QA{xi)h{x 1 ) + h(x' 1 )Q A (x 1 ) N(N-l) ^^ 

V. A = T A n = N 1 VA,a{Xl,X2)g(Xi,X2) (31) 



rl A = — [h A (xi) + h^x-y) J H gA,A\x\,x<i) H 9A,a{xi,x 2 ) . (32) 

a^A 

The expectation value of energy, Trl/jH^}, using Eq. (32) for a time-independent density, reduces to the time- 
independent result [Eq. 26]. 

Note that T-L A is a non-Hermitian operator with real eigenvalues. That the eigenvalues of the operator are real is 
demonstrated by the fact that the expectation value of H A is real for any state of the system. The non-Hermitian 
character of the T-L A results from the fact that the subsystem, as an open system, evolves nonunitarily in time [14]. 
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IV. APPLICATION TO ENERGY TRANSFER 



The above results are applicable to molecular systems of any constituency, and with components at any intermolec- 
ular distance. Here we apply the definition of the subsystem Hamiltonian to donor-acceptor systems and illustrate the 
agreement between results calculated based on Ha and well known limits that are applicable when the subsystems 
are suitably separated. 

Energy Transfer problems are solved in different limiting regimes that are classified according to the spatial proximity 
of, and coupling strength between, interacting molecules. We will follow the development of these limits as outlined 
by Parson [15]. To illustrate analytically the agreement of our result with the Forster and Dexter limits, consider 
a sample calculation in the limit of weak coupling and insignificant subsystem overlap. The composite system is 
partitioned into two distinct parts, A and B with normalized subsystem eigenstates |aj) and \bj). To simplify the 
algebra we assume that A and B are far enough apart that the overlap, (a,i\bj), is sufficiently small so as not to 
require rcnormalization of the product wavefunction without necessarily being zero. We denote Hamiltonians of the 
two isolated systems by Ha and He] by definition then (ai(l)|-ffyt(l)|oj(l)) = E ai 8ij, and similarly for B. 



A. Two-electron Case 



For A and B sufficiently separated the system Hamiltonian of the non-interacting composite system is given by 
H = Ha + Hb and the eigenstates of the composite system are \ipij) = -^\ai(l)bj(2) — bj(\)cii(2)) . A general state, 
"f, of the system is then a superposition = J^ij c ul 7 / , u)> with the density matrix given by p = \^){^\. Using the 
subsystem Hamiltonian "Ha the energy of subsystem A, Ea, in state "J is 



E A = (*|Wa|*> 

= ( Yj C ^ 



/ ijkl 



(33) 



In the limit of well separated subsystems the wavefunction of the system vanishes at the boundary between the 
subsystems. As discussed in the preceding section on time evolving states we are therefore justified in using Ha as 
given in Eq. (26). Focusing on the (ij,kl) matrix element gives 



{ipij\H A \^ki) = A^i|0A(l)Ml)|^w> 



N(N — 1) 



(^■|^e A , Q (i,2) 3 (i,2)|^) 



(34) 



Moreover 



(iMMi) = = ri<Oi(l)&i(2) - & i (l)o i (2)|ex(l)^(l) 



1 

71 



^(1)6,(2)1^(1) (35) 



and 



(^■1^(1,2) = (^-16^(1,2)5(1,2) = i(a i (l)6 J -(2)-6 J -(l)a i (2)|e A (l)6 A (%(l,2)=0 (36) 



(1>ij\9A,B{l, 2) = (^|e^ B (l,2) 5 (l,2) 

4(a i (l)&i(2)-6 i (l)o i (2)|i( 

' (^(1)^(2) - 6,(1)^(2)1.9(1,2) (37) 



= -^(a i (l)6 i (2)-6 i (l)o i (2)|-(e A (l)e B (2) + e B (l)e A (2)) ff (l J 2) 



2V2 



Hence, 



(A 3 \U A \ipki) = ^(1)6,(2) - 6 i (l)a i (2)|0 A (l)W|o fc (l)6,(2) - 6i(l)o fc (2)) 

= y(a 4 (l)6 J -(2) - & 3 -(l)a ( (2)|e >1 (l)ft(l)|a fc (l)& J (2) - bi{l)a k (2)) 



N(N 



-^5>,(1)&,(2) 



2 

a 

- 6 i (l)a j (2)|e XiB (l,2) 5 (l J 2)|a fc (l)6,(2)-6i(l)ofc(2)) . (38) 
With the number of electrons N = 2, and the summation over partitions A and B, 



(iPt 3 \H A \ipki) = (a i (l)6 i (2)|ft(l)|a fc (l)6 I (2)-6i(l)o fc (2)> 

1 X 1 

+ 2^2 



a i (l)6 i (2) - 6 J -(l)a i (2)| fl -(l,2)|o fc (l)6z(2) - 6j(l)o*(2)) 



= fa I ft ( 1 ) I a fc ) + \ {Oj bj I g{ 1 , 2) I a fc 6; ) - § (a^- 1 g ( 1 2) | k a k ) (39) 

EA t Sik Jijki Kijki 

where the Coulomb interaction matrix element, Jijki , can be (i) the matrix element, Jtjij, for interaction between 
charges at A and B, (ii) the matrix element, Jijik(Jijkj), for interaction between a charge at A(B) with a transition 
density at B(A), and (iii) the matrix element, Juki, f° r interaction between transition densities at A and B. The 
coupling, (iii) above, between the transition densities at A and B is dominant in the Forster limit and it mediates 
the motion of Frenkel excitons between subsystems. When the separation of subsystems A and B is large relative 
to the spatial extent of either subsystem the transition density coupling can be calculated accurately in the dipole- 
dipolc limit [ ]. The Exchange interaction matrix element, K^ki, accounts for the interactions between densities, 
Pij = cubj, at A and B. The overlap between eti and bj decays exponentially with distance from the boundary of each 
subsystem. Therefore, the Exchange interaction can dominate only if the subsystems are very close to one another or 
if the singlet-singlet Coulomb interaction is symmetry forbidden. This constitutes the Dexter limit of energy transfer. 
Thus 

E A = Tr{pH A } = c* 1 c k i(ipij\T A 'H\'4>ki) 

ijkl 

— \ 2_j c ij c kl(EAi$ij,kl + Jij.kl — Kij,kl) 
ijkl 

= taf^ + 5 ]C c *i c ki{ J ij,kl - K ijM ) (40) 

ij ijkl 

Hence, the energy of subsystem A is a sum of the energy of the isolated molecule plus half of the coupling energy, 
comprising the coulombic coupling J, and the exchange coupling K . Since the interaction energy is split equally 
between the two subsystems the total electronic coupling energy of the system is twice the interaction energy assigned 
to subsystem A. The electronic coupling matrix element is the sum of the Coulombic coupling, J, and the Exchange 
coupling, K, 

ycoul = j + R _ ( 41 ) 

These are the Forster and Dexter coupling results obtained by Parson [15]. 

B. Multi-electron Case 

To extend the analysis of the previous section to multi-electronic systems, let n(M) be a single determinant M- 
electron wavefunction of subsystem A. 

|H(M)) = -J=A\a?(l)c%(2) ■ ■ ■ o° (M)) , (42) 

where the anti-symmetrization operator 

A = £(-l) p P (43) 
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sums over all possible permutations P of the M electrons 1, 2, . . . , M. The factor (— l) p represents the parity of the 
permutation P. Similarly T(K) is the antisymmetric K-electron wavefunction of subsystem B. 

|r(K)) = -^=A\b\ (1)&£(2) ■ • • b T K {K)) . (44) 

The antisymmetrized N-electron wavefunction, ^(N), of the composite system is then 

l*(N)> = ( KKN ^ K)l ) 1/2 F\U(M))\T(K)) , (45) 



where F sums over the permutations between the first M and the remaining K electrons. Then 



l*<N)', = \l Kl(N m K)l F^=A\aT(l)a$(2) ■ ■ ■ a*(M))-?=A\l$(l)$(2) ■ ■ ■ b T K {K)) 

^F{A\a?(l)aK(2) ■ ■ ■ a* (M ))}{A\% (1)^(2) • • • b T K {K))} (46) 

In accord with our definition, the energy of subsystem A is given by 
E A = Tr{pH A } 

= jV(*|Ml)W + N(N ~ 1) (*lgA,A(l,2)lg) + N(N ~ 1) -£(*\g Ata (l,2)\V) 

M 



^(^(1)^(1)1^(1)) 
i=l 
M 

+ 2'<o i (l)a i (2)| fl (l,2)|a i (l)o i (2)> - (a j (l)o J -(2)| fl (l,2)|o i (l)a < (2)) 



i=l 
j>i 

M K 

+ J2J2« 1 ^)\9(h2)\a i (l)b j {2)) - (a,(l)6 J -(2)| 5 (l,2)|6 i (l)a i (2)) 

= T A + V a ,a + \v a ,b (47) 

In the limit where A and B are infinitely separated Va,b = so that limfl^oo E A = T A + Va,a we properly recover 
the energy of the isolated subsystem A. 



V. SAMPLE COMPUTATIONAL RESULTS 



We present numerical results for coherent energy and charge transfer dynamics in several elementary systems: (a) 
the hydrogen molecule, H2, where each atom is considered as an open subsystem, and (b) two hydrogen molecules 
interacting with one another. Issues of nuclear antisymmetrization are neglected. In the former case the interacting 
subsystems have open shell electronic configurations while in the latter case the subsystems have closed shells. Hartree- 
Fock (HF) molecular orbitals are obtained in both cases using Gaussian 03 [17] at various bond lengths R, for which 
EET results are shown. When the molecule is in a stationary vibrational state, EET results would be a weighted 
average over those shown below for various R values. Despite the simplicity of the systems, the results are enlightening. 



A. Interacting Hydrogen Atoms 

Figure 1 presents the coherent time evolution of energy and population between two interacting hydrogen atoms 
initially in the state \^sa2pz,b) & t various fixed internuclear separations. In each panel of the figure the curves near 
the top show the time evolution of the population on each hydrogen atom, together with their time average. The 
curves near the bottom of each panel show the time evolution of the energy of each atom, together with their time 
average. 
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a) 



b) 



H Bond Length, R = 1.0 a.u. 



H 2 Bond Length, R= 1.4 a.u. 
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SO 60 70 

Time (a.u.) 



c) 



d) 




H Bond Length, R- 10.0 a.u. 



- -0.2 a 



40 50 60 
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FIG. 1. H2 molecule prepared in electronic state |1sa2p z s). Coherent energy and charge transfer dynamics as function of time 
for bond lengths (a) R = l.Oa.u., (b) R — lAa.u., (c) R = l.Sa.u., and (d) R — lO.Oa.w. 



The figures confirm several intuitive expectations. At large separation (Fig. Id), when the atoms are nearly isolated, 
there is negligible exchange of both energy and charge. One hydrogen atom has the electronic ground state energy of 
Isa, while the other hydrogen atom has the electronic excited state energy of 2ps- As the bond length shortens, with 
an accompanied increase in the electronic coupling strength, we observe the expected increase in the rate of energy 
transfer (Fig. l(a-c)). 

The figure also reveals two interesting features. First, energy transfer is here seen to be perfectly correlated with 
charge transfer for all R. Also, the amount of energy and charge exchanged coherently between the hydrogen atoms 
maximizes at some intermediate separation, consistent with the fact that no energy or charge will be exchanged in 
the two limits, i.e., in the limit of the unified atom and that of infinitely separated atoms. 

To further illustrate these patterns, consider next a hydrogen molecule prepared in a superposition of its ground 
and first excited electronic state (Fig. 2). The dominant contribution to the ground electronic state is from the <j\ s 
configuration. The first excited electronic state has a primarily ai s a* s character, where <j\ s ~ (ls^ + lsg) is the sigma 
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a) b) 



H Bond Length, R = 1 .0 a.u. H Bond Length, R = 1.4 a.u. 




10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 

Time (a.u.) Time (a.u.) 

FIG. 2. H2 molecule prepared in electronic superposition state |*) = l/\/2(|^ 'ground) +|(/>i»* excit ed)) ~ l/v / 2(|o"i s ai s ) + |ai s CTj' s )). 
Coherent energy and charge transfer dynamics as function of time for bond lengths (a) R = 1.0a. u., (b) R = lAa.u., (c) 
R = l.Sa.u., and (d) R = 10.0a.it. 



bonding molecular orbital while a* s ~ (Is a — Isb) is the sigma anti-bonding molecular orbital. A superposition of 
these states gives rise to coherent dynamics between the Is orbitals on the two hydrogen atoms. Here too, the rates 
of charge and energy transfer as the two hydrogen atoms approach each other increase, and there is an optimal R at 
which energy and charge transferred are maximized. At infinite separation limit the two hydrogen atoms have equal 
energy, corresponding to that of isolated |ls) hydrogen atoms. 

The dynamics in this case, which involves the superposition of the two lowest singlet electronic states look much 
simpler than that in Fig. 1. By contrast, the state in Fig. 1 projects onto several CI eigenstates of the H2 molecule 
and the dynamics therefore involves multiple timescales corresponding to many eigenenergy differences. 
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FIG. 3. A sketch of the ground and excited electronic potential energy surfaces of the H2 molecule 



Decoherence 



As a further example we would like to consider the effect of decoherence, e.g., decoherence induced by collisions 
of gaseous H 2 molecules, on the coherent energy and charge transfer dynamics. We assume that the collisions as 
elastic, resulting in a dephasing of otherwise coherent dynamics. Moreover, we write the complete wavefunction of 
our system as a product of the electronic and nuclear wavefunctions in the Born-Oppenheimer approximation. For 
the superposition state studied above (Fig. 2), we may write the wavefunction as: 

\9(r,R)) = -^=\* g (r,R)) + ^=\* e (r,R)) 

= ^=\%(r;R))\ Xg (R)) + -^\Mr;R))\Xe(R)) (48) 

where \Xg(R)} an d \Xe(R)} are taken to be the ground vibrational wavefunctions on the respective ground and excited 
electronic potential energy surfaces (Fig 3). The electronic energy on hydrogen atom A is then: 

(E A ) = (*(r,R)\H A \V(r,R)) 

= \ Yi I dRxt(R)Xi(R) ( f drift (r;R)H A 1>i(r;R) 



i=g,e 



J dR Xg {R)Xe{R) (J 



drip*(r; R)R A ip e (r; R) 



(49) 



In order to make the expression more concise we define the following quantities 

LUi — Ei + (n, + ^)hvi 

uj eg = (E e - E g ) + (n e + -)hv e - (n g + \^)hv g 



Em = 



1 



dR X *(R)x i (R) / dr^(r;R)H A Ur;R) 



A A = 3tj dR X * g (R)Xe(R) (J 



dnl>*Jr;R)H A Mr;R) 



(50) 



where i = g,e. Moreover, the time dependence of the wavefunction is easily included; only the cross term will be 
time dependent. Thus, the time dependent energy on Hydrogen atom A in the concise form 

(E A (t)) = 1/2 (E Ag + E Ae ) + A A cos(uj eg t) 
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FIG. 4. Decoherence induced dynamics in (a) Energy, and (b) Charge of the the H2 molecule prepared in a superposition of 
the two lowest electronic states 



The net result is a time-dependent site electronic energy, reflecting the time-dependent superposition state. 

In addition, in the presence of a collisional environment the vibronic energy levels of the system are expected to 
fluctuate about their mean positions leading to different Hydrogen atoms in an ensemble of molecules acquiring an 
arbitrary phase. The expectation value of the time-dependent energy on hydrogen atom A can then be expressed as 



(E A (t)) = § (E Ag + E Ae ) + A A J dwcos{wt)f(u) 

Here f(co) describes the distribution in energies due to environmental collisions. We assume a normal distribution of 
phases with mean uj eg . In our simulations the standard deviation is taken to be 10~ 6 Hartree (corresponding to a 
collision time on the order of 25 ps). The results of the simulation are shown in Figure 4. We observe that for the 
given parameters the electronic energy transfer dynamics fades out on a timescale of 4 x 10 6 a.u. (100 ps). 



B. Interacting Hydrogen Molecules 

As a second example, consider the energy and population transfer dynamics between two interacting Hydrogen 
molecules at various separations. As an example, we fix the bond distance within each hydrogen molecule at the 
equilibrium bond length of 1.4a. it. and vary the distance between the centers of the two H2 molecules. For each 
nuclear geometry the Hartree-Fock molecular orbitals of the system are determined using the 6-31G basis set, and 
the singlet excited states of the system are determined at the Cl-Singles (CIS) level. 

Figure 5 presents the coherent time evolution of energy and population between the interacting molecules when the 
system is initiated in a superposition of its ground and first excited singlet CIS electronic state. Similarly, Figure 6 
shows the dynamics for a superposition of the ground and a higher lying electronic excited state. The curves for the 
evolution of population and energy are respectively near the top and bottom of each panel. 

The dynamics of the interacting closed shell subsystems can be compared to the dynamics of the interacting open 
shell system presented earlier. As before, the amplitude of population and energy exchanged between the interacting 
subsystems decays with increasing separation. The correlation between the population and energy dynamics within 
each subsystem, and the anticorrelation between the dynamics of the two subsystems also persists. However, two 
important differences can be discerned. First, note that charge transfer is only significant at separations below 5.0 
a.u., where as significant energy transfer may persist beyond 10 a.u. This is in contrast to the two hydrogen atoms 
where we found that charge transfer always accompanied energy transfer. Moreover, because we no longer have the 
unified atom limit, the amplitude of energy and population transfer no longer shows a clear trend at small separations. 



14 
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b) 



H2 — H2 separation, R = 2.0 a.u. 




H2 — H2 separation, R - 3.0 a.u. 



I -2.4 
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H2 — H2 separation, R = 5.0 a.u 
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H2 — H2 separation, R = 10.0 a.u 
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FIG. 5. Interacting H2 molecules prepared in an electronic superposition state of the ground and first excited CIS state. 
^} = 1/V2( I <t>ground) + \4'i ! ' t excited})- Coherent energy and population transfer dynamics as function of time for bond lengths 
(a) R = 2.0a.u., (b) R = S.Oa.u., (c) R = 5.0a.u., and (d) R = lO.Oa.u. 



VI. SUMMARY 



We have introduced a well defined operator Ha that allows, given a time dependent or time independent wave- 
function |\&(t)), the computation of electronic energy on a local site in a composite molecular system, as Ea — 
(^/(t)\'H.A\^(t)). The definition resolves numerous problems arising from electron interchange antisymmetrization 
that exist with naive approaches. Further, the computation of Ea, being an integral over an operator that is a func- 
tion of the electronic degrees of freedom, automatically includes decoherence effects due to other degrees of freedom 
in the molecule. The resultant operator has been shown to give appropriate results in various limits and to provide 
insight into electronic energy dynamics in small molecular systems. Applications to larger systems arc underway. 
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FIG. 6. Interacting H2 molecules prepared in an electronic superposition state of the ground and third excited CIS state. 
^} = l/\/2( 1 4>ground) + \<Pyd exC i te d) ) ■ Coherent energy and population transfer dynamics as function of time for bond lengths 
(a) R = 2.0a.u., (b) R = S.Oa.u., (c) R = 5.0a.u., and (d) R = lO.Oa.u. 
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